Manuscript submitted to 
AIMS' Journals 

Volume X, Number OX, XX 200X 



Website: http://AIMsciences.org 
pp. X— XX 



INTERACTIONS OF POINT VORTICES IN THE 
ZABUSKY-MCWILLIAMS MODEL WITH A BACKGROUND 

FLOW 



COLM CONNAUGHTON 

Warwick Mathematics Institute and Warwick Centre for Complexity Science 
University of Warwick 
Coventry, CV4 7AL, United Kingdom 

John R. Ockendon 

Oxford Centre for Industrial and Applied Mathematics 
Mathematical Institute, University of Oxford, 
24-29 St Giles', Oxford OX1 3LB, United Kingdom 



Abstract. We combine a simple quasi-geostrophic flow model with the Zabusky- 
McWilliams theory of atmospheric vortex dynamics to address a hurricane- 
tracking problem of interest to the insurance industry. This enables us to 
make predictions about the "follow-my-leader" phenomenon. 

1. Introduction. 

1.1. Insurance industry motivation for the "follow-my-leader" problem. 

The 2005 Atlantic hurricane season is famous as the most active and expensive 
Atlantic hurricane season since records began. In addition to thousands of deaths, 
damage to property and infrastructure was estimated to have amounted to 130 
billion USD. The Mexican states of Quintana Roo and Yucatan and the U.S. states 
of Florida and Louisiana were each struck twice by large hurricanes. These included 
Katrina, the most expensive natural disaster in the history of the United States. 
The insurance industry is naturally interested in the question of whether this is 
mere coincidence or whether there are correlations between the tracks of intense 
hurricanes. Risk estimation models used in the industry often treat the probability 
of hurricanes making landfall in a particular area as independent Poisson processes 
characterised by their historical mean. The possibility of even weak correlations, 
particularly between large storms, may be an important source of systematic error 
in these models. 

This question was posed by Lloyds and explored during the 73rd European Study 
Group with Industry which took place at Warwick in April 2010 [5]. Statistical 
analysis of the historical data was done which did not prove conclusive. In tandem 
with this statistical analysis a more basic fluid dynamics question was studied: how 
do the pairs of vortices embedded in a larger scale "steering flow" influence each 
other when their separations are large? 
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Our work will build on the basic theories of oceanic-scale atmospheric flows 
described in [3, 7] and the theories of hurricane dynamics proposed in [1, 2, 8]. We 
will rely heavily on the atmospheric vortex dynamics model proposed in [9] and the 
observational studies reported in [4]. 

1.2. Inviscid vortex dynamics in 2-D: summary of the results of the Study 
Group. The most elementary theory of 2-D vortex dynamics concerns the flow that 
results when vortices of strengths I\ at Xj(£) = (xi(t), j/t(t)), i = 1 . . . n, move in a 
background potential flow, 

dip dipo" 



1 dx ' dy J \ dy ' dx 

where ipo(x, y, t) is a prescribed streamfunction. The theory asserts that the stream- 
function is 

■4>{x,y,t) = ip (x,y,t) + V — ^ logr-i, (1) 

* — ' l IT 

1=1 

where rf = (x — Xi) 2 + (y — y,) 2 , and that the (xi(t) , yi(t)) evolve according to the 
Hamiltonian system: 



dxi dip (x.) 



dt dy 
dyi 5-0o (x) 



dt dx 



x= Xj (t) J=1 27Tr ij 



r,-(t) (»(*) -%•(*)) 

(2) 



9-Trr 2 

x=x,(t) J = 1 Z7rr ij 



Eq. (1) ensures that the flow is potential flow away from the vortices and Eq. (2) is 
the Hclmholtz condition that each vortex moves with the velocity that would have 
existed in its absence. This condition can be justified by smearing the vorticity into 
small patches (see [8] for the use of this idea to model hurricanes) around (xi,yi) 
and considering the momentum balance for these patches, which reveals that no 
relative velocity can exist between the patch and its ambient free stream. 

A very crude model for the "follow-my-leader" problem is to consider the motion 
of two equal vortices placed in a steering flow comprising of a uniform potential 
flow impinging upon a solid wall (see Fig. 1(A)). The flow in the upper right 
quadrant may be thought of as representing the Atlantic anticyclone and the wall 
as representing a blocking pattern in the Gulf of Mexico. We took 

*l>o = Uxy, (3) 

The streamlines of the large scale steering flow are plotted in Fig. 1. The only 
interesting feature is that the flow has a hyperbolic point at (0, 0). With r = the 
point vortices simply follow the streamlines shown in Fig. 1(A). 

We performed a number of numerical experiments to show how r ^ affects this 
passive advection of vortices. Fig. 1(B) shows the tracks of two equal vortices having 
r = 1 which start from the same point as in Fig. 1(A), (this was (4.00, 0.01) in the 
figure) but separated by a delay in time. The lead vortex follows the same path until 
the second vortex is introduced. The lead vortex then undergoes a deflection which, 
because it occurs when the lead vortex is near to the stagnation point, subsequently 
results in a completely different southwards trajectory to the northwards trajectory 
it would have followed in the absence of the second vortex. Note that the second 
vortex also experiences a deflection but still follows roughly the same path. 
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(A) Unperturbed track (B) Track perturbed by a second vortex 
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Figure 1. Motion of point vortices in the steering now, Eq. (3). 
(A) shows a single vortex following the streamlines of the steering 
flow. (B) shows the large change in trajectory which can occur if 
the original vortex (solid circles) is near a stagnation point when a 
second vortex (open circles) is introduced at a later time. 

Of course this large deviation in the trajectory of the first vortex does not occur 
for all configurations of the two vortices. The conclusion to be drawn from this 
simple model is, as one might expect, that the question of whether vortices follow 
each other is not straightforward even in this simplest case. Trajectories which 
come close to stagnation points of the steering flow are very difficult to predict. 



1.3. Rotational background flow. A serious complication arises when we try 
to generalise the above approach to rotational steering flows, U = (Uq, Vq) where 
V • Uo = and (U • V) U = — Vp in dimensionless variables. Now Uo is derived 
from a streamfunction that satisfies 

J[^ o ,A^o]=0, (4) 

where J [/, g] is the Jacobian operator, 

j \f g ] = - ^1 2l 

' dx dy dy dx 

Considering just one moving vortex for simplicity, it is now easy to see that if we 
suppose that 

Y 

tp(x, y, t) = ip Q {x, y, t) + j— log [(a; - x ) 2 + (y - y ) 2 ] , (5) 

then 

Aip = Aip + T S(x - x ) S(y - y ). 

Also 

dip dtpo r y-yo 



dy dy 2 ir (x - x a ) 2 + (y - y ) 2 

dip dipo T x — xo 

dx dy 2ir (x - x ) 2 + (y - y a ) 2 
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and 



dx 
dAiP 

dy 
dAip 
dt 



A 



A 



dx 
dy 



TS'(x - x ) S(y - 
TS(x - x )S'(y 



Vo) 

yo) 



Tx 6'(x - x ) 6{y - y ) + Fy S{x - x ) 5'{y - y )- 



Hence when we collect terms in the vorticity equation, 

dAiP 



dt 



+ J[Vo,A^ ] =0, 



(0) 



we find that the coefficients of S'(x — xq) S(y — yo) and S(x — xo),5'(y ~ yo) are 
Xq + and yo — g^f respectively, in accordance with the Helmholtz condition, 
Eq. (2). However there are also terms like 

(X - Xq) 6(X - Xq) 5' (y - t/p) 

(x - xq) 2 + (y - y ) 2 
which vanish and then terms like 

y - yo d x - x d 



(x - x ) 2 + (y- y ) 2 dx (x - x ) 2 + (y - y ) 2 dy 
which only vanish in a potential flow. What has happened is that the introduction 
of the vortex has interfered with the global distribution of vorticity meaning that 
we have to solve the full vortical Eulcr equation, Eq. (6) everywhere. Since almost 
no steering flow for hurricane tracks can be approximated by a potential flow, this 
means that the results of point vortex models such as that presented above have 
very limited usefulness unless the generation of background vorticity caused by the 
introduction of the point vortex is small in some sense. This may be the case if the 
strengths of the moving vortices are weak so that in Eq. (5) can be replaced by 
a small parameter, e. If we try the ansatz, 

<P(x, y, t) = Mx, y,t) + e [log [(x - x ) 2 + (y - y ) 2 ] + + 0(e 2 ), (7) 

where ipi is non-singular at (xo,yo), then we find that we can make the terms of 
O(e) in Eq. (6) if tpi satisfies: 

dtpi 



Of 



ch/jo A d_ _ 
dy dx 

y-yo 



di/jp A d 

dx dy 
d 



dipi d dipi d 
dy dx dx dy 

X — Xq d 



_(x - x Q ) 2 + (y ~ y ) 2 dx (x - x ) 2 + (y - y ) 2 dy 

so that the size of the perturbation can, in principle, be obtained by solving a linear 
equation. 

The feedback from the vortex dynamics to the steering flow is even more dramatic 
when a vortex of strength T is introduced at the origin into a weak steering flow 
with streamfunction e^Jo(x 1 y). Then a simple perturbation argument shows that 
the streamfunction, eijj(,x,y,t) evolves according to 



A 



dip _ r 

~dt ~ 2ir(x 2 +y 2 ) 



d_ 

dx 



d_ 
dy 



At/) 



on a timescale which is short compared to that of the vortex motion. Hence the 
vorticity in the steering flow has to be rearranged before the vortex trajectory can 
be computed. 
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This discussion illustrates one of the main difficulties of building a more realistic 
model of the steering of hurricanes by large scale atmospheric motions since such a 
model must necessarily incorporate the rotation of the Earth and the background 
vorticity which it induces. The problems associated with the incorporation of the 
Earth's rotation can be partially addressed using a model originally devised by 
Zabusky and McWilliams [9] although we note that we are still required to assume 
that the production of background vorticity remains small. We discuss this model 
in the next section. 



2. Zabusky-McWilliams model of point vortices on the /3-plane. 

2.1. The Charney equation. The atmosphere is, to leading order, in geostrophic 
balance. That is to say, the largest terms in the equations of motion which approxi- 
mately balance each other, are the horizontal pressure gradient and the coriolis force 
induced by the Earth's rotation (see [7] for a full discussion) . Any reduced model of 
atmospheric dynamics must take rotation into account if it is to have any chance of 
being applicable. The most basic pde model of the quasi-two-dimensional dynamics 
characteristic of large scale atmospheric motion are the so-called quasi-geostrophic 
equations (again discussed in great detail in [7]), which, in the simplest case of 
pure barotropic motion can be reduced to the Charney equation [3] for a single 
scalar streamfunction, ^(x, t), of a two-dimensional spatial coordinate, x = (x,y) 
and time, t. Written on the so-called /3-plane, where x denotes the longitudinal 
direction and y denotes the latitudinal direction, the Charney equation is written 

r) 8ty 

-(Ay- 1 2 y)+p— + J[M>,AM>}=Q, (8) 

where 7 is the inverse of the Rossby deformation length, and /3 is the rate 
of variation of the Coriolis parameter with latitude. The geostrophic velocity is 
obtained from the streamfunction by taking the curl of <j)z: 

The Charney equation only differs from the two-dimensional Euler equations by the 
addition of the two linear terms. The physics is somewhat different however. In 
particular, in addition to playing the role of the streamfunction for the geostrophic 
velocity, $ has a direct physical meaning. Because the atmosphere is assumed to 
be shallow, it is proportional to the hydrostatic pressure and, equivalently, to the 
deviation of the depth of the atmosphere from its equilibrium depth. Some simple 
manipulations allow it to be written as a Lagrangian conservation law: 

— [A*- 7 2 *-/3t/l = ^ =0. (10) 
Dt L J Dt y ' 

The quantity, Q, which is conserved along fluid trajectories is called the potential 

vorticity. The j3 and 7 2 terms add two important ingredients to the basic dynamics 

described by the Euler equations. The j3 in the potential vorticity, means that 

the intrinsic vorticity of a fluid parcel changes when it moves in the y (latitudinal) 

direction. This can be shown to produce a restoring force on fluid displacements 

in the latitudinal direction which introduces waves into the model. These waves, 

known as Rossby waves are an important feature of large scale atmospheric motions. 

The 7 2 term introduces a potential energy penalty for the generation of large values 

of or, equivalently, for large deviations of the atmospheric thickness from its 
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equilibrium depth. This term sets the observed characteristic scale for large scale 
motions in the atmosphere. 

2.2. Zabusky-McWilliams model. The idea of Zabusky and McWilliams [9] was 
to try to find the analogue for the Charney equation of the point vortex represen- 
tation of the Euler equations. The issue is not straightforward since the /3 term is 
intrinsically continuous and does not lend itself easily to a discrete representation. 
The basic idea of Zabusky and McWilliams was to start from the vortex circulation: 

fi(x,i) = A*-7 2 *, (11) 

and discretise it on a set of N time-dependent vortices, {xj(i) = (#,(<), Ui(t)), i = 1, . . . , N}, 
as one would do for the Euler case: 

N 

A* (x, f) - 7 2 *(x, t) = ^( x - ( 12 ) 

i=l 

As discussed in section 1, the points Xj move with the fluid motion: 

cfajj _ 9*(x,t) 
(if 9y 

cfy, _ 0¥(x,t) 



(13) 



x=x t (t) 



eft <9x 



=Xi(i) 

The important difference, however, is that the vorticity is not conserved along 
fluid trajectories on the /3-plane. Rather the potential vorticity is conserved as in 
Eq. (10): 

Ki(t)+0yi(t) = const = ql 0) (14) 

where is the potential vorticity of vortex i at t = 0, which reveals the depen- 
dence of the vortex strength on position. The streamfunction is obtained from the 
potential vorticity by inverting the modified Hehnholtz operator in Eq. (11) and 
using Eq. (12) and Eq. (14): 

*(x,t) = -^Jdyn(y,t)K Q (<y\\x-y\\) 

1 N 

= -^Y,^ ~ PVi(t)) K ( 7 \\ X - Xi (t)\\), (15) 

Z7T * — ' 

t=l 

where Kq(z) is the Bessel function of the second kind of order zero. Eq. (15) and 
Eq. (13) now yield closed equations for the time evolution of the vortex centres. 
The conservation of potential vorticity enormously changes the dynamics compared 
to the case of Euler point vortices [1, 8], leading to quite complicated trajectories 
even in the case of equal strength vortices. 

2.3. Zabusky-McWilliams model with a background flow. We now discuss 
one way of adding a background flow, ^o( x )! to the Zabusky-McWilliams model. 
We represent the contribution from a discrete set of point vortices by ip(x,t), so 
that: 

*(x,t) = ¥ (x)+V(x,t). (16) 
Likewise the potential vorticity can be decomposed into an ambient part, Qq(x), 
and a discrete part, q(x,t): 

Q(x,t) = Q (x) + flf(x,t), (17) 
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where 



)(x) 



N 



g(x,t) = A^(x,i)- 7 2 </;(x,i) =5j/s i (t)*(x-x i («) 



(18) 
(19) 



Noting that for atmospheric flows, the total potential vorticity, Q(x, £), should be 
conserved following each vortex: 

Q(x i (t),t) = Q(x i (0),0). (20) 

Integrating this equation over the interior of an infinitessimal contour enclosing x$ 
we obtain 

AQ (xi(t)) + **(*) = AQ ( Xi (0)) + Ki(0) (21) 

where ^4 is the area enclosed by the infinitessimal contour. From this, we obtain 
an approximate evolution equation for the vortex intensities in the presence of the 
background flow: 



Ki(t) = Ki(fl) + A [Q (xi(0)) - Q (xi(t))] . 



(22) 



In Sec. 1 we have noted the complication caused by the coupling between the 
vortex dynamics and the evolution of the vorticity in the steering flow. Here we 
will circumvent this difficulty by assuming that the circulations associated with 
the vortices are small compared to the circulation in the initial steering flow as in 
Eq. (7). Hence, to lowest order, the motion of the vortex centres are again obtained 
from Eq. (13): 



~dt 

dyi 
dt 



(x) 



dy 

dx 



0V(x,t) 



x=x»(*) d v 

dx 



(23) 



x=Xj(i) 



X=Xi(t) 



x=x 4 (t) 

It remains to express the discrete part of the streamfunction, ip(x,t) in terms of the 
positions of the vortex centres. This is done, as before, from Eq. (19). We obtain 



1 N 

-Y 

2tt ^ 



Ki(t) X ( 7 ||x-X l (0||). 



(24) 



with Ki(t) given by Eq. (22). The equations of motion resulting from Eqs. (23) and 
(24) written out explicitly are: 



dxi 
lit 

dyi 

dt 



a* (x) 



dy 
5* (x) 



dx 



x= x< (t) £j 2 7 r|x i (t)-x i (i)| 

V «i(t) ^l(7 \Mt) - Xi(t)l) &•(*) - Z,(i)) 



x=x;(i) 



E 

i=i 



2^ | Xi (*)-x^)| 



-, (26) 



where the Ki(t) are expressed in terms of the x,(f) through Eq. (22). These equations 
differ from the original Zabusky-McWilliam model in two respects. Firstly the 
equations of motion contain a velocity coming from the background flow so that a 
single point vortex will move following the streamlines of the background flow even 
though that flow is vortical. Secondly, the modulation of the vortex intensities is 
now depends on the background flow as specified by Eq. (22) . This ensures that the 
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model remains consistent with the principle of conservation of potential vorticity. 
Clearly these equations reduce to the original model if the background flow is absent. 

We have not yet discussed possible forms for the background flow, \E'o(x). We 
require it to be a stationary solution of Eq. (8) or, equivalently, Eq. (10). There 
are a large number of stationary solutions of the Charney equation. When ^ is 
independent of time, Eq. (10) gives 

j[*,A*- 7 2 *-/3y] =0. (27) 

From this, is is clear that 

A* - 7 2 * - /3 y =F(V) (28) 

yields a solution for any function F(*5/). In this paper, we are mostly interested in 
the interaction between point vortices in the presence of a background flow rather 
than in the details of the background flow itself. For this reason, we consider here 
two specific forms of the background flow: 

1. Uniform zonal current 



This corresponds to a background flow consisting of a uniform westerly 
flow. Streamlines are straight lines. The corresponding streamfunction is 

*o(x,y) = -Uy. (29) 
2. Inertial boundary current 



An inertial boundary current [7] occurs when a uniform westerly flow from 
x = oo encounters a straight north-south boundary at x = 0. It is the analogue 
for Eq. (8) of the well-known potential solution of the two dimensional Euler 
equation describing a uniform flow impinging upon a flat plate. It has a 
stagnation point at the origin. The streamlines are exponential curves and 
the streamfunction is 



$> (x,y) = Uy 



1 — exp 




(30) 



3. Numerical results. 



3.1. Nondimensional equations. Let us measure lenghts in units of the Rossby 
deformation length, 1/7 and velocities in terms of the characteristic velocity, U 
of the background flow. The natural unit of time is then (<yU)~ . Introducing 
dimensionless variables, x', if and defined by 

x=-x\ t= — t', and* = — 

7 7 U 7 

the nondimensional version of Eq. (8) is (we immediately drop the primes): 

^(A*-*) + ^ + J[*,A9] =0, (31) 

where 

A 

7 2 t/ 

is the dimensionless /3-parameter. Our model background flows are 
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Physical quantity 


Value 


Notes 


Rossby deformation length, 7 _1 


10 e m 


1000 km 


Beta parameter, /3 


1.6 x KT 11 m-^-s- 1 




Hurricane force wind velocity 


33 ms~ 1 


118 km/hr 


Typical radial extent of hurricane force winds 


1.6 x 10 5 to 


160 km 


Typical hurricane eye radius 


2.4 x 10 4 to 


24 km 


Typical hurricane core area, A 


1.81 x 10 9 to 2 




Typical circulation of a hurricane 


3.34 x 10 7 m 2 s- 1 




Horizontal scale of steering flow 


5 x 10 6 to 


5000 km 


Typical velocity of steering flow 


8ms _1 


17 mph 


Traversal time 


6.26 x 10 5 s 


s» 7 days 



Table 1 . Physical values of parameters relevant to hurricane dy- 
namics taken from [6] . Hurricane circulation has been estimated as 
r = 2 tt Rv based on hurricane force winds of velocity v — 33tos _1 
extending to a distance R = 1.6 x 10 5 to. The vortex core area, A, 
is estimated based on the eye radius. 



1 — exp 



fix 



(Uniform zonal current) (33) 
(Boundary current) (34) 



The dimensionless versions of Eqs. (25) and (26) are 



dxi 
~dt 



9*o(x) 



dy 



N 



dt 



d^o (x) 
dx 



X=X;(t) j = l 

N 



y, R 3 {t) K^jt) -x 3 -(t)|)( yi (f) - Vj {t)) (; , V! 



2tt \ Xi (t)- Xj {t)\ 



^ 27T|x i (t)-X J -(t)| • 1 ' 



X=Xj(t) 

Here the dimensionless circulations of the point vortices are 

Ri(t) = Ri(p) + A [Q (xi(0)) - Qo(xi(t))] , (37) 

with A a dimensionless area (which can be absorbed into the background flow) and 
Ki(0) is the dimensionless initial strength of vortex i: 

«i(0)7 



Ki(0) = 



U 



Qo is obtained from Eq. (33) or Eq. (34): 

Q (x) =AV -y +h- 



(38) 



(39) 



The upshot of all of this is that the only control parameters in the problem are 
the initial intensities, Ri(0), of the point vortices relative to the strength of the 
background flow. In what follows we shall take all vortices to have equal initial 
strength. 

Some geophysically plausible values for the various parameters in the original 
equations are presented in table 1. The corresponding dimensionless values used in 
the numerics are summarised in table 2. 
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Dimcnsionlcss parameter 


Value 


Dimcnsionless beta parameter, /3 
Initial vortex circulation Ri(0) 
Dimcnsionless vortex core area, A 
Dimensionless traversal time 


2.0 
4.2 
2.0 x 10~ 3 
5.0 



Table 2 . Values of dimensionless parameters used in the numeri- 
cal simulations obtained from the physical values presented in table 
1. Note that the vortex circulation, although large compared to 
U /7 is less than that of a typical Atlantic anticyclone circulation 
(~ 10 to s ), in accordance with the comments made before 
Eq. (23). 




2 4 6 8 10 12 

lime delay (days) 



FIGURE 2. Deflection after 11 days of a vortex moving in a back- 
ground current caused by a second vortex originating at the same 
location a specified time later. Open circles correspond to a uni- 
form backround current, solid circles to a separating background 
flow. 

3.2. Deflection of one vortex by another with background current. Fig. 2 

shows the results of a numerical experiment designed to quantify the effects of 
vortices on each other. The parameters used are as in Table 1 and the axes have 
been re-expressed in physical units for ease of interpretation. We used two different 
background flows, a uniform current flowing to the left (west) (open circles) and 
separating background flow corresponding to an inertial boundary current (solid 
circles). A single point vortex simply follows the streamlines of the background 
flow as it would have done in the classical theory of potential flows. In order to 
quantify the influence of other vortices on this trajectory we started a vortex off at 
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(A) Unperturbed track (B) Delay of 0.5 days 




-1 012345 -1 012343 



(C) Delay of 1 day (D) Delay of 2 days 




-1 012345 -1 012345 



(E) Delay of 3 days (F) Delay of 4 days 




-1 012345 -1 012345 



FIGURE 3. Tracks of vortices in a separating background flow for 
different delay times. Panel (A) shows the unperturbed track in 
the absence of a second vortex. Very large deflections, of the order 
of the system size, are observed for delays of up to 4 days. 



a chosen point in space and then inserted a second vortex at the same point after a 
given time delay of r. We then compared the position of the first vortex at a later 
time, T = 11 days, with the position it would have been at had the second vortex 
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not been inserted. The follow-my-leader effect can then be quantified by plotting 
this deflection of the first vortex from its unperturbed path as a function of the 
time delay between the vortices. Clearly, for large delays, the deflection should be 
smaller as the time delay increases. 

For the case of a uniform background flow, it decreases exponentially as the 
delay increases. We see that for a delay of 7 days, typical for Atlantic hurricanes, 
the deflection is about 10 km. Since this is less than the core radius of our vortices, 
below which our model cannot be considered meaningful, this means that there is 
effectively no interaction between the vortices for such separations. Of course, much 
stronger deflections of the order of hundreds of km can be observed for vortices 
separated by 2 days or less which is reflecting the well-known fact that vortices 
interact very strongly when they get close to each other. 

For the separating background flow, the typical deflection is orders of magnitude 
larger even though all model parameters remain the same. The reason is clear from 
Fig. 3 which illustrate the tracks followed by the two vortices for different values 
of the delay. One sees that the second vortex can easily deflect the first one onto a 
track which subsequently diverges from the initial track due to the presence of the 
stagnation point in the flow. 

4. Conclusion. We have used a quasi-geostrophic model together with a physically 
plausible law of vortex dynamics to model hurricane tracks in the presence of an 
ocean-scale steering flow. Our model has enabled us to quantify the sensitivity 
of the "follow-my-leader" phenomenon to the presence of stagnation points in the 
steering flow. It would be of interest in the future to compare the predictions of 
our simulation with the probabilistic models that are used in risk estimation. 
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